Mixed state non-Abelian holonomy for subsystems 
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Non-Abelian holonomy in dynamical systems may arise in adiabatic transport of energetically 
degenerate sets of states. We examine such a holonomy structure for mixtures of energetically 
degenerate quantal states. We demonstrate that this structure has a natural interpretation in terms 
of the standard Wilczek-Zee holonomy associated with a certain class of Hamiltonians that couple 
the system to an ancilla. The mixed state holonomy is analyzed for holonomic quantum computation 
using ion traps. 
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I. INTRODUCTION 

Non-Abelian holonomy in simple dynamical systems 
has found application in a variety of contexts, such as 
nuclear rotations of diatoms 0, 01 j molecular Kramers 
doublets i i, nuclear quadTu^ole resonance (NQR) 
0, 0, II] , dynamics of deformable bodies 0, , the n- 
body problem ^3, semiconductor heterostructures |12| |. 
and ion traps Recently, it has been pointed out 

fl^ that non-Abelian holonomy may be used in the 
construction of universal sets of quantum gates for the 
purpose to achieve fault tolerant quantum computation. 
This has triggered further work on holonomy effects for 
quantum computation [HMHE^j^^ i21j, IH HI III 
and quantum information |25l I26ll27l l28l| . 

An important issue for holonomic quantum informa- 
tion processing is its robustness to imperfections, such as 
decoherence and random unitary perturbations. The ef- 
fect of imperfections on the quantum gates has been ana- 
lyzed 17,, .29, 30, 31j, supporting the alleged robustness of 
holonomic quantum computation. It remains however to 
address the appearance of non-Abelian holonomy related 
to the motion of the quantal states themselves. Here, we 
wish to address this issue by putting forward an attempt 
to define a concept of non-Abelian holonomy in terms of 
incoherent mixtures of energetically degenerate quantal 
states. 

The main idea of this paper is to relate the mixed state 
time evolution of a system S to the holonomic transfor- 
mation of a certain class of pure state systems consisting 
of S coupled to an ancilla system A. For the whole sys- 
tem S + A, this coupling creates in a purely geometric 
fashion an entangled state whose partial trace over A cor- 
responds in general to a nonpure state of S. From this 
perspective, we may regard the present analysis as the 
identification of a class of S — A interaction Hamiltonians 
having degenerate eigenspaces that admit a non-Abelian 
holonomy structure for the mixed states of S. This struc- 
ture reduces to the Wilczek-Zee pure state holonomy jl| 
of the iS* system when the S~ A interaction vanishes. We 



examine the present type of mixed state holonomy in the 
case where the pure state limit corresponds to the ion 
trap system in Refs. ,13. .IGJ . 



II. MIXED STATE HOLONOMY 

Let a quantal system S be exposed to the Hamiltonian 
H(q), q being some external control parameters that vary 
around a closed path C : t € [0,T] —f qt in parameter 
space. Let Hs be the Hilbert space of S and let M{q) 
denote the linear span of the iV-fold degenerate instanta- 
neous eigenstates l^ail)), a = 0, . . . , N — 1, N < dim 7^5, 
of H{q). Under the condition that T is large enough 
so that the adiabatic theorem holds, any initial state 
|'0(O)) e M{qo) is mapped to G M{qt). Thus, 



l^(o)> = \U<io)) 

\-<P{t)) ^ exp ( -i I E{qt)dt 



(1) 



with U{t) unitary, E(qt) the common instantaneous en- 
ergy of the degenerate subspace, and we have put h = 1 
from now on. It follows that 



U[C] = Pexp (^-^A 



(2) 



is the holonomy transformation for C. Here, P is path 
ordering and A = ^fj^Ai^{q)dq^ is the connection one- 
form with 



Aab,M = {Uq)\d^Uq)) 



(3) 
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the components of the Wilczek-Zee gauge potential and 
= d/dq^,. 

We now apply the above Wilczek-Zee prescription in 
order to achieve the main idea of this paper, which is to 
propose a concept of mixed state non-Abelian holonomy 
for subsystems. To this end, we first identify a proper set 
of degenerate entangled energy eigenstates by consider- 
ing the extended Hilbert space H. = Hs ® Ha , where Ha 
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is an ancilla Hilbert space taken to be identical to Tis and 
the vectors in 7i being purifications of the mixed states 
of S. The essential criterion is that this set should de- 
fine the standard pure state holonomy of the S system in 
the product state limit. This may be achieved by letting 
H{q) be the A^-fold degenerate Hamiltonian of the com- 
posite system S + A and let the members of the instan- 
taneous set of orthonormal eigenvectors span- 
ning the degenerate subspace of H{q), be iso-entangled, 
i.e., that they correspond to purifications of mixed states 
Tr^|Sa((7))(Sa(g)| all having the same eigenvalues A^, 
fc = 0, . . . , iV - 1. If A; 1 and Xk^i 0, then it fol- 
lows that {\Ea{q))} -> {\^a{q))<»\f)}, for some \ip) € Ha, 
which establishes the desired connection to the holonomy 
in the pure state case. 

Now, the iso-entangled energy eigenstates may explic- 
itly be taken as 

AT-l 

|Sa(g)) = \/^\^la+kUi<l)) \Vk), (4) 

k=0 

where [a + k]N = a + k mod N. For this orthonormal set, 
we obtain the gauge potential as 

Aab,p. = (St(g)|(9^|Sa(g)) 

= J2 ^''^^l'>+k],^{q)\d^\C[a+k]^iq))- (5) 

k=0 

The resulting unitarity U[C] for the loop C in parameter 
space is the desired mixed state non-Abelian holonomy 
of the set of energetically degenerate reduced density op- 
erators Tryi|Sa((z))(Sa(g)| pertaining to the subsystem S. 

To further analyze the physical meaning of the pro- 
posed form of mixed state holonomy, let us note that the 
Hamiltonian corresponding to this setup may be written 
as 

H{q) = Eiq)J2\E,iq)){Eaiq)\+AH{q) 

a 

EE HE{q) + AH{q), (6) 

where E{q) is the energy of the degenerate subspace and 
AH{q) acts only on the orthogonal complement. Thus, 
H{q) induces interaction between S and A along the path 
C. For example, in the TV = 2 case, we have 

He = £;(g)(a^(Q)® [Ao|(po)(</?o| +Ai|^i)(^i|] 

+ y/%>^af{q)^at), (7) 

where 

^o^(?) = Mq)){ioiq)\ + Mq)){U<i)l 

^Uq) = Mci)){(i{q)\ + Mq))Mq)\, 

= \ipo){ipi\ + \ipi){ipol (8) 

In the pure state limit AqAi ^ we have \Ea{q)) = 
1^0(9)) <8) I'Po), Va, which decouples the Hamiltonian and 



leads to the pure state Wilzcek-Zee gauge potential in 
Eq. @. 

Now, since S in general interacts with A, one should 
expect that the evolution of the density operator p is 
nonunitary. To check this, one may assume that p{0) = 
TiAHmWO)], with IV'(O)) = \Ea{qo)) e ns®-HA. 
calculate the final density operator as 

p{T) = TvaH{T)){^{T)\ 

= 5^C/:,[C]C/„,[C]Tr^|S,(0))(Sb(0)|, (9) 

b,c 

and compare the initial and final degree of mixing for in- 
stance by computing Trp^. Due to the appearance of 
off-diagonal terms in p(T), it is indeed apparent that 
Trp2(r) ^ Trp2(0), in general. 

III. PHYSICAL EXAMPLE 

A scheme to realize geometric quantum computation 
for ion traps has been proposed [T ^ (se e also Ref. 
Here we extend the U2 gate in Ref. [l^ , using our gauge 
potential for mixed states. 

Suppose we have an ion with three metastable ground 
states denoted by |0), and \a), and one excited state 
|e). A laser couples the ground states with the excited 
state, as described by the Hamiltonian 

H = \e){{0\LJo + (l|c^i + {a\uja) + h.c, (10) 

where ujQ,uji,uja are complex coupling parameters that 
may be varied by slowly switching the laser field on 
and off. To realize the generalization of the U2 gate 
of Ref. restrict to paths that correspond to the 

parametrization lvq = lu sin 6 cos (f>, ui — lu sin 9 sin and 
Wa = w cos 6. Under this condition, there is a pair of zero 
energy dark states of the form 

\DQ{6,(j))) = cos6'cos0|O) -f cos0sin0|l) — sin0|a), 
\Di{0,(j>)) = -sin(/)|0) -HCOS0I1). (11) 

For simplicity, in the adiabatic limit loT 3> 1, we may 
assume that the system's density operator diagonalizes 
in this dark state basis, i.e., 

p{0,<l>) = ^{l + r)\Do{9,cl,)){Do{e,c^)\ 

+ i(l-r)|i?i(0,0))(Z?i(0,<^)| (12) 

with — 1 < r < 1. By using Eq. jsj, we obtain 

A = iray cos 6d<j), (13) 

where a^Tf^y, <^z are the standard Pauli operators in the 
|0), |1) basis. However, this connection one-form is only 
valid for a patch excluding the north and south pole of the 
sphere, where the coordinates are not uniquely defined. 
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This can be seen by considering an infinitesimal loop 5C 
around the north pole, say, yielding 



U[5C] = Pexp 



A 



sc 



-iray27T 



To get around this problem we make the gauge transfor- 
mation 



A 



Am = VmAV, 



' N 



with 



Explicitly, 



V, 



N 



corresponding to a product state in Hs <S) Ha- Con- 
versely, one may transform a pure input state into a 
mixed one using the mixed state holonomy U[C]. To 
(14) illustrate this latter point, let {ip'^) = + r)(|0) + 

^ l^o) + ^\/(l^(|0) - ^ and m = 

iy(rT7)(|o)-|i))®i(po)+iy(r^(|o)-f|i))^i^i), 

so that p'{0) = |0)(0| is obtained from the purification 
[■^(O)) = --i=(|7/ig) + I'ip'i)). After having traversed the 
loop C, we obtain 

mr)) = 4^[(c/oo[C] + t/lo[C])|v^) 



A 



N 



-iray{\ 



(15) 
(16) 
(17) 



which is a valid gauge potential for the northern hemi- 
sphere. 

Now, since [AN-tpiO)^ AM:,j,{0')] — for any pair 9,6' 
along C, path ordering is not needed and the mixed state 
holonomy becomes 



U[C] 



(18) 



where is the solid angle enclosed by the path C in 
parameter space. In the pure state limit r =^ 1, U[C] 
reduces to the U2 gate in Ref. [l6l| . 

As mentioned above, we expect the mixing of the state 
to change around the loop C. To check this, let p(0) = 
i(l -hr)|0)(0| + 5(1 - r)|l)(l|, with purification |V'(0)) = 

^|(1 + 0|0) ® |<^o) + ^Kl-r)|l) ® 1^1) = liJo). After 
having traversed the loop, we obtain 



|V(r)> = ;7oo[C]|Vo) + c/oi[c]|^i) 

— cos{ril)\'i(jo) + sin{ril)\tlJi) 



(19) 



with IVi) = + <E> 1^0) + ^i(l-r)lO) <S> 

which is a purification of 



p{T) = ^(^I + sin{2rn)ax+rcos{2rn)a, 



(20) 



whose Bloch vector traces out an ellipse in the {x, z) plane 
with half axes (1, |r|) when Q, varies. We obtain 

Trp2(T) = i(l + r2 + (l-r2)sin2(2rr!)), (21) 



which equals Trp2(0) = i(l 



if and only if |r| = 1 



(pure states) or = n7r/(2r), n integer. 

we obtain that the initial mixed 
state p(0) polarized along the z axis has transformed into 
the pure state p(r) = \{I±a^). This may be understood 
by looking at the final purification, which for n = reads 



V2 



V2 



(|0) + |1)) 



+ ^\{^-r)W2)) (22) 



V2 

+ (C/oi[C] + C/n[C])|^l) 

cos(rrj + 7r/4)|V'o) 
-Kcos(rr2 - 7r/4)|V''i), 



(23) 



which is a purification of the final density operator 

p'{T) -r sin{2rn)<7^ + cofi{2rn)a^^ (24) 

that is nonpure for \r\ < 1 unless fl — nn/ {2r), n integer. 
Note that the accessible final states again constitute an 
ellipse in the {x, z) plane, now with half axes (|r|, 1). 

More general manifolds of final states are obtained if 
the initial S + A state is mixed. As an illustration of this, 
consider the input state 

^^(0) = ^mii'ol + ^^l^iX^il (25) 

with —1 < i? < 1 and \ipo)A^i) defined above. This 
state evolves into g{T) = U[C]g{0)W[C], yielding the fi- 
nal density matrix of S as 

p{T) = i + i? sin (2rf7) + rR cos (2rf7) cr^) , (26) 

which corresponds to a shrinking of the ellipse in the 
{x, z) plane by a factor \R\. 



IV. CONCLUSIONS 

Holonomy effects for mixed states, as first conceived by 
Uhlmann js^l and later reconsidered in Ref. [s^l , have re- 
cently attracted interest due to its potential importance 
to fault tolerant quantum information processing. These 
previous analyzes have mainly been concerned with the 
Abelian case, while, on the other hand, universal sets 
of quantum gates must involve noncommuting one- and 
two-qubit gates. In other words, to examine the effect of 
nonpure input states in all-geometric quantum informa- 
tion processing requires a concept of non- Abelian mixed 
state holonomy. 

In this paper, we have introduced such a concept of 
non-Abelian holonomy for adiabatic transport of ener- 
getically degenerate mixed quantal states. The holon- 
omy effect may be understood as arising in certain types 
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of slowly changing Hamiltonians that couple the consid- 
ered system to an ancilla. 

The present type of holonomy effect has an inherent 
additional richness compared to the standard Wilczek- 
Zee pure state holonomy in that it may change the mixed- 
ness of the input state. We have explicitly demonstrated 
that this feature may transform a mixed state into a pure 
one and vice versa by purely geometric means. This pro- 
vides a method for robust coherent control of quantal 
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